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ABSTRACT 

In this paper, we recompute contributions to the spectrum of the nonlinear integrated Sachs- 
Wolfe (iSW)/Rees-Sciama effect in a dark energy cosmology. Focusing on the moderate non- 
linear regime, all dynamical fields involved are derived from the density contrast in Eulerian 
perturbation theory. Shape and amplitude of the resulting angular power spectrum are similar 
to that derived in previous work. With our purely analytical approach we identify two dis- 
tinct contributions to the signal of the nonlinear iSW-effect: the change of the gravitational 
self-energy density of the large scale structure with (conformal) time and gravitational lenses 
moving with the large scale matter stream. In the latter we recover the Birkinshaw-Gull effect. 
As the nonlinear iS W-effect itself is inherently hard to detect, observational discrimination be- 
tween its individual contributions is almost excluded. Our analysis, however, yields valuable 
insights into the theory of the nonlinear iSW-effect as a post-Newtonian relativistic effect on 
propagating photons. 
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1 INTRODUCTION 

The tiny fluctuations in the temperature of the cosmic microwave 
background (CMB), its primary anisotropics, reflect the state of 
the Universe approximately 400,000 years after the big bang. Pho- 
tons emanating from the last scattering surface, however, interact 
on their way to today's observer with the intervening large scale 
structure (LSS), thereby altering this picture of the early Universe. 
One of these secondary anisotropics is the integrated Sachs-Wolfe 
(iSW) effect (Sachs & Wolfe 1967). The interaction of CMB pho- 
tons with time-evolving gravitational potentials in the LSS causes 
additional fluctuations in the CMB temperature. While being blue- 
shifted when entering a potential well the photons also experience 
a redshift on their way out of the potential. A net change in temper- 
ature results when blueshift and redshift do not compensate. 

Being a secondary effect, the signal of the iSW-effect is quite 
small in comparison to the CMB temperature itself. Its importance, 
however, arises from its dedicated sensitivity to cosmic fluids with 
non- vanishing equation of state (Crittenden & Turok 1996) pro- 
moting it as a valuable tool in the investigation of dark energy and 
non-standard cosmologies (Lue et al. 2004; Zhang 2006), as well as 
the Universe's spatial curvature (Kamionkowski & Spergel 1994). 

Exploiting its strong correlation with the galaxy distribution 
the iSW-effect has been detected with high significance by various 
groups (Boughn & Crittenden 2004; Vielva et al. 2006; McEwen 
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et al. 2007; Giannantonio et al. 2008). The constraints they derive 
for the matter density Cl m and the dark energy equation of state 
W give support to the ACDM model, which describes a spatially 
flat universe mainly constituted by cold dark matter (CDM) and 
cosmological constant A. 

The nonlinear iSW- or Rees-Sciama (RS-) effect (Rees & 
Sciama 1968; Seljak 1996; Cooray 2002; Schafer & Bartelmann 
2006) refers to any temperature anisotropy caused by nonlinearly 
evolving gravitational potentials. In contrast to the linear iSW- 
effect it is a small scale phenomenon dominating the total iSW sig- 
nal on multipoles I > 100 but is nevertheless hard to detect because 
of the primary CMB (Cooray 2002; Schafer et al. 201 1). 

Resulting from nonlinear structure formation the statistics of 
the RS-efFect is necessarily non-Gaussian. The non-Gaussianities 
induced by the RS-effect, however, are small. Spergel & Goldberg 
(1999) showed that the CMB bispectrum due to the RS-effect is 
undetectable. Also mixed polyspectra of the nonlinear iSW-effect 
temperature perturbation and the tracer galaxy density field, as re- 
cently investigated by Jurgens & Schafer (2012) in continuation of 
the work of Schafer (2008), do not reach sufficiently high signal- 
to-noise ratios. Therefore, stacking methods (Granett et al. 2008) as 
an alternative to a statistical detection of the RS-effect are of great 
interest. 

Seljak (1996) and Cooray (2002) provided two different per- 
turbative descriptions of the nonlinear iSW-effect yielding compa- 
rable results. In this work we take up the ansatz of Cooray (2002). 
While his results make extensive use of findings from the halo 
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model, we are more interested in an analytical treatment. There- 
fore, we concentrate our analysis on the weakly nonlinear regime 
where the statistics of the underlying fields are not too far from be- 
ing Gaussian and second order perturbation theory is well applica- 
ble. We compute the nonlinear iSW-effect using solely the density 
field and its derivatives. Although the applicability of this ansatz is 
naturally limited, it allows for deeper insights into the underlying 
physical processes. Our approach is well suited for identifying the 
distinct contributions to the nonlinear iSW-effect as well as eluci- 
dating their physical origin. 

This article is structured in the following way: We begin with 
a short compilation of the key formulae describing growth and 
statistics of the density contrast in dark energy cosmologies (Sec- 
tion 2). In Section 3 we compute the spectrum of the nonlinear 
iSW-effect. Furthermore, we derive the resulting two-point statis- 
tics and present the corresponding angular power spectrum. We ex- 
tend our analysis in Section 4, where we discriminate and interpret 
the different contributions to the nonlinear iSW-effect. In Section 5 
the physical interpretation is continued and deepened by pointing 
out analogies with classical field theory as well as the theory of 
gravitomagnetic potentials. Finally, we summarize our results in 
Section 6. 

Throughout this work we choose a spatially flat vvCDM cos- 
mology as reference. The initial perturbations in the CDM com- 
ponent are assumed to be adiabatic and Gaussian distributed with 
variance <x 8 = 0.8, where we set the spectral index n s to unity. The 
matter content is described by D ra = 0.25 and Q.\, = 0.04 while the 
value of the Hubble constant is set to H = 100/jkm/s/Mpc with 
h = 0.72. 



2 COSMOLOGY 

2.1 Dark energy cosmologies 

In a spatially flat Friedmann-Lemaitre-Robertson- Walker universe 
the time evolution of the scale factor a is governed by the the Hub- 
ble function H(t) = dloga/d?: 



H 2 (a) 
"1 



■ Sl m a 3 + (1 - Q m )exp 



Jd log a 
a 



[1 + w(a)] 



(1) 



The expansion of the homogeneous background is thus completely 
described by the matter content fl m (in units of the critical density) 
and the equation of state of the dark energy fluid w(a). Following 
Chevallier & Polarski (2001) the latter can be phenomenologically 
parametrized as 



w(a) = wo + (1 - a)w a , 



(2) 



covering a huge variety of different dark energy models. Setting 
Wo = - 1 and w a = one recovers the equation of state of a cosmo- 
logical constant. 

Since photons follow null geodesies, conformal time dq = 
dt/a and comoving distance x can be used interchangeably allow- 
ing for a simple relation to the scale factor: 



A ■ 



la * 



da 



H(a) 



(3) 



Naturally, (comoving) distances are measured in units of the Hub- 
ble distance xh = c /Hq which sets the scale up to which Newtonian 
gravity is applicable. 



2.2 CDM power spectrum and linear structure growth 

The fluctuations in the cold dark matter (CDM) component are de 
scribed by the density contrast 



Six) 



with 



Pa 



^mPcrit- 



(4) 



In linear theory <5 is a statistically homogeneous and isotropic Gaus- 
sian random field and thus completely characterized by its power 
spectrum 

(6(k)S'(k')) = (2n) 3 S D (k - k')P S6 (k). (5) 

For the power spectrum we choose a power law modulated by an 
appropriate transfer function, i.e. 



for which we use the fit proposed by Bardeen et al. (1986) 



(6) 



T(q) 



log(l + 2.34g) 



2.34a 

+(6.71 ? ) 4 ) 4 . 

In this formula the wavenumber q 
parameter 

T = £l m h exp 



(l + 3.89a + (16.1o) 2 + (5.46a) 3 



(7) 



n h l + 



V2fr 



kT is rescaled by the shape 



(8) 



(9) 



in order to account for a non-vanishing baryon density Db 
(Sugiyama 1995). Finally, the amplitude of the power spectrum is 
set by the variance of the density contrast <x 2 smoothed on the scale 
R = 8Mpc//i 

al = J ^k 2 P 6S (k)W 2 R (k) = dlogk A 2 (k)W 2 R (kR), 

where W R (y) = 3j\ (>')/(>') is the Fourier-transform of a spherical 
top hat function {ji(x) denoting the f-th spherical Bessel function 
of the first kind, Abramowitz & Stegun 1972). 

During linear evolution the time dependence of the density 
contrast is completely encapsulated in the growth function D + (q) so 
that S(k, rj) = D + (q)6o(k) (normalized to unity today). It is obtained 
by solving the growth equation (Turner & White 1997; Wang & 
Steinhardt 1998; Under & Jenkins 2003) 

d 2 1 / dlog#\ d 3 

— D + (a) + - 3 + — 2— —D + {a) = —Q. m (a)D + {a). (10) 
aa 1 a \ dloga / da 2a 2 



3 NONLINEAR INTEGRATED SACHS-WOLFE EFFECT 

The relative change in temperature CMB photons experience when 
they traverse time-evolving gravitational potentials on their way 
from the last scattering surface to today's observer along the direc- 
tion h is given by the line of sight integral (Sachs & Wolfe 1967) 



AT iS w(") 2 r XH d 
J cmb c Jo drj 



(11) 



The gravitational potential is directly related to the density fluctua- 
tions via the (comoving) Poisson equation 



8(x,Tj) 



AO(jc,i/) = -H 2 £l, 

2 a 



(12) 



Its time derivative can be most easily accessed in Fourier space 
where the Poisson equation is readily inverted 
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dr. 



d 5(k,t]) 



k 2 dr, 



(13) 



In the derivation of the linear iSW-effect one usually inserts 
6(k,a) = D+(a)So(k) revealing that during matter domination 
where D + (a) = a the iSW-effect vanishes. For our purposes, how- 
ever, it proves advantageous to carry out the time derivative for- 
mally 



a or, 



(14) 



Obviously, the time evolution of the potential is sourced by two 
different effects: on the one hand the expansion of the homogeneous 
background, represented by the first term, and on the other hand by 
the time varying structure growth. As mentioned before, in epochs 
dominated by a fluid with vanishing equation of state (w = 0) these 
two effects just balance. 

Since we are interested in nonlinear contributions to the iSW- 
effect we express the time derivative of the density fluctuations by 
means of the continuity equation. 



■^-6(x,r]) + di\j(x,rj) ■ 
or, 



0. 



(15) 



Due to matter conservation, the change of the dark matter density 
in time needs to be compensated by the divergence of the corre- 
sponding flux 



j(x, r,) = [1 + 6(x, ?])]v(x, rj). 



(16) 



In the linear regime 6 <k 1, thus the momentum density 
6(x,r,)v(x,r,) is negligible and one can immediately read off the 
familiar Fourier representation of the peculiar velocity field 

.dD + {rj) k 



v(k, r]) 



J d 3 xv(x,r])e' k 



dr, k 2 



So(k). 



(17) 



In Fourier space, the dark matter current density becomes a convo- 
lution 



j(k, if) = v(k, r,) + 



C d 3 k' 
J (2^r? 



v(k- k',r,)6(k',r,), 



(18) 



so that the time evolution of the gravitational potential finally can 
be written as a function of density contrast and peculiar velocities 

^-®(k,ri) 
dr. 



^H 2 n m ^[H(r])6(k, rj)- ik- (v(k, 77) 
+ / §^ v(k - k '^ 6(k '^)l 



(19) 



Due to the symmetry of the convolution the expression for the cur- 
rent density is not unique. For example, one could have explicitly 
symmetrized equation (18). However, for our purpose, aiming at 
the different contributions to the RS-effect, the form given above is 
most useful as will become clear below. 

So far, our derivation is essentially the same as the one pre- 
sented by Cooray (2002). Cooray (2002), however, focused on the 
deeply nonlinear regime and proceeded by using approximations 
resulting from the halo model approach. In contrast to this, we are 
interested in an analytical description of the contributions to the 
nonlinear iSW-effect which arise from dark matter currents. There- 
fore, we restrict our analysis to the moderate nonlinear regime. Ob- 
viously, in this regime the simple relation between peculiar veloci- 
ties and density contrast, given in equation (17), is no longer strictly 
valid. The CDM particles cease to exclusively follow the gradients 
of the potential. However, in order to get an analytical estimate of 
the nonlinear contributions, it is still a useful approximation. Thus, 
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Figure 1. Growth function of the linear (dashed line) and nonlinear (solid 
line) iSW-effect along with linear structure growth (dotted line), for varying 
equation of state parameters of the dark energy fluid. 



we take nonlinear corrections to the continuity equation, manifest 
in the momentum density, into account, but we maintain the linear 
properties of the underlying fields, namely density contrast and pe- 
culiar velocities. Consequently, by virtue of equation (17), there is 
only one dynamical field, namely the density contrast, involved in 
our analysis. Equation (18) then may be written as 

dD + {rj) f d 3 k' k-k! 



j{k,ri) = -iD + (rf)- 



th] 



I 



(2n) 3 \k - k'\ 2 



S (k - k')S (k') (20) 



omitting its linear part. This expression for the current density 
has already been successfully used in the computation of the 
power spectrum of the Ostriker-Vishniac effect (Ostriker & Vish- 
niac 1986; Jaffe & Kamionkowski 1998; Castro 2003). Since in our 
framework all nonlinear contributions to the iSW-effect are solely 
sourced by the convolution part of the flux it is convenient to define 
its rescaled divergence 

-PA 



0(M) = - 



D + (rf) dD + (jj) 



r d 3 k' 

J On) 3 



k(k 



dr] J (2;r) 3 k 2 \k- k 
whereas the linear contributions are captured by 
6 D + {rf) S (k) 



k -l6 (k-k')6 (k')(2\) 



m,v) = 



dr, 



(22) 



With this the change of the gravitational potential takes the compact 
form 



^-O(A:,77) = -^ 2 fi 1I 
or, 2 



[f(k,r,)-&(k,v)]. 



(23) 



In Q m /xli we recognize the gravitational coupling constant medi- 
ating between the evolving matter fields and the resulting change 
in the gravitational potential. Equations (21) and (22) underline the 
different time evolution of the linear and nonlinear iS W-effect. For 
illustration we plot the relevant combinations of the growth func- 
tion and its derivative for three different equation of state parame- 
ters of the dark energy component in Figure 1 . 

In order to characterize the statistics of the source field 
r(k, r,) = j^Q>(k, r,) we compute its power spectrum next 



{T(k)r'(k')) = (2n) 3 S D (k - k')P TT (k). 



(24) 



Because we assume the density contrast being a Gaussian random 
field its bispectrum vanishes and the power spectrum completely 
separates in its linear and first order nonlinear contributions 
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Prr(k) = P„(k) + Pee(k). 



(25) 



The linear contribution is given by Pff(k) = P M (k)/k , while for the 
evaluation of the second term we invoke Wick's theorem yielding 



Pee(k) 



r d 3 k' 
J 



k(k 



k 2 \k-k'\ 2 
xP 6S {\k-k'\)P 65 (k'). 



k') k ■ k' 
+ 



k 2 k'- 



k (k- k') 
k 2 \k-k'\ 2 



(26) 



At this point we would like to emphasize that our signal of the non- 
linear iSW-effect completely originates from the trispectrum of the 
density field. In contrast to this, Cooray (2002) entirely neglects 
this term using arguments derived under the halo approach. Fur- 
thermore, he does not take terms involving cross-power spectra be- 
tween the density and velocity fields into account. In our approach, 
however, these terms do contribute since they are contained in the 
trispectrum. Another difference with respect to the work of Cooray 
(2002) results from its somewhat artificial discrimination between 
velocity-density correlations on the one hand and correlations of 
the density field and its (conformal) time derivative on the other 
hand. For the latter, Cooray (2002) borrows once more a result from 
the halo model assuming that for the fields involved the Cauchy- 
Schwarz inequality actually becomes an equality. In our approach, 
however, this far-reaching assumption can be avoided. 

The computation of the convolution integral in equation (26) is 
most readily carried out in a coordinate system where the wavevec- 
tor k is aligned with the x-axis and by introducing spherical coor- 
dinates such that k! = ak(ji, sin r? cos 0, sin i? sin 0) with a = k' Ik 
and f-i = cos &. Expressing the CDM power spectrum in its dimen- 
sionless form, we find 



\ r°dtt r 1 

2 Jo 



{aii-\)[a{2n 2 -\)-n] 



2 Jo a 2 J_, r ^(a 2 - 2afi + 1) 7/2 
xA 2 (k Vo- 2 - lap + l) A 2 (ak), 



(27) 



where we have denned Ag (fc) = k 3 P @ e(k)/2n 2 . The corresponding 
angular power spectrum is obtained by Limber projection (Limber 
1953; Bartelmann & Schneider 2001) 



n r x 
e Jo 



X d X W 2 xx (x)A 2 x > 



(k = e/x) . 



In case of the linear iSW, the weight function is given by 
d D + (a) 



w„(x) = 



dn 



(28) 



(29) 



(30) 



while for the first order nonlinear one we have 
a arj 

In Figure 2 we plot both weighting functions for varying equation 
of state parameter of the dark energy fluid. As expected, the weight 
function of the nonlinear iSW-effect peaks at comoving distances 
where the one of the linear iSW-effect has already started to de- 
cline rapidly. At high redshifts Cl m (a) approaches unity so that here 
growth function and scale factor coincide. Consequently, the linear 
iSW-effect vanishes identically (cf. equation 22). 

The resulting power spectra are shown in Figure 3. Reflecting 
the additional fluctuation amplitude of the matter distribution on 
small spatial scales, the spectrum of the nonlinear iSW-effect is 
much flatter than the linear one. In case of the ACDM model the 
nonlinear signal surpasses its linear counterpart at angular scales 
smaller than about one degree, i.e. C > 100. For time evolving dark 
energy models this crossing has shifted to smaller scales. However, 
being of second order in the density contrast, the amplitude of the 
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Figure 2. Weighting function of the linear (dashed line ) and nonlinear (solid 
line) iSW-effect for three different values of the equation of state parameter 
of the dark energy component. 
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Figure 3. Angular power spectra of the linear (dashed line) and nonlinear 
(solid line) iSW-effect for three different values of the equation of state 
parameter of the dark energy fluid. 



nonlinear spectrum is much less sensitive to the underlying dark 
energy equation of state than the linear one. Remarkably, the spec- 
trum of the nonlinear iSW-effect presented here is very similar to 
the result of Cooray (2002) regarding its shape and amplitude, even 
though Cooray (2002) worked in the deeply nonlinear regime while 
our analytical estimate is dedicated to the translinear regime. 

It is interesting to also investigate the mass dependence of 
the nonlinear iSW-effect. We content ourself to the ACDM model 
and smooth the power spectrum of the density contrast on different 
mass scales by introducing a Gaussian filter function 

P S6 (k) — » PsAk)S 2 R (k) with S R (k) = exp(-k 2 R 2 /2). (31) 

The smoothing scale R is set by the corresponding mass scale 
Smooth = ^PaitQmP 3 ■ In Figure 4 we show the angular power 
spectra of the linear and nonlinear iSW-effect for three different 
masses, namely 10 11 , 10 12 and 10 13 solar masses. Excluding ob- 
jects with masses larger than 10 13 solar masses the power decreases 
considerably on small scales (C > 1000), demonstrating that mostly 
cluster-size objects contribute to the signal. 
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Figure 4. Contributions of different mass scales to the angular power spec- 
tra of the linear (dashed line) and nonlinear (solid line) iSW-effect. 



4 CONTRIBUTIONS TO THE NONLINEAR 
INTEGRATED SACHS- WOLFE EFFECT 

A deeper analysis of equation (21) reveals some physical insights 
into the origin of the nonlinear iSW-effect. Mainly, one can identify 
two distinct sources. First, we rewrite equation (21) as the sum of 
two new fields 



9(4,17) i 
with 

and 

0*g(M) 



® SG (k,ri) + ® BG (k,ri), 



D + (U) dD + (l) 
a dr. 



f d'kf S ( 
J (2tt) 3 \k 



(k - k) 



k'\ 2 



So(k') 



(32) 



(33) 



D + (rj) 3D + (ij) 



f d 3 k' k ■ k' 
J (2^F \k-k'[ 



■6 (k - k')So(k'). (34) 



dr, J (2tt) 3 \k - k'\ 2 

Concentrating on the first term, we proceed by making use of the 
Poisson equation and obtain 



©sc(M) 



H 2 ^ 



a dr. 



J (^ Oo 



(k - k')S (k') 



— +aH(a)\ 
dr] 



likewise in position space 
8 



®sg(x,t}) 



+ aH(a) 



<b(x, 7,)S(x, J]). 



(35) 



(36) 



Keeping in mind that we have already absorbed the inverse of 
the prefactor in the definition of the weighting function (cf. equa- 
tion 30), the product <b(x)6(x) can be interpreted as the gravita- 
tional self-energy density associated with the density contrast and 
its own gravitational potential. Consequently, the field ® S a may 
be interpreted as the (conformal) change of the gravitational self- 
energy density of the LSS together with the variation arising from 
the homogeneous expansion of the background. 

The interpretation of the second term is more involved. First, 
we rewrite equation (34) in a more suggestive form 



®BG(k,ri) 



r d 3 k' 

J (2n) 3 ' 



, D + {rf) S (k - k') dD + {r,) k 



a \k- k'\ 2 



dr, k 2 



So(k'). (37) 



Using the Poisson equation, we recognize that the first term has 
the form of a gradient of the potential. The second term, however, 
resembles the velocity field (cf. equation 17). Thus, we are tempted 
to write 



® BG (x, rf) =: -2V0)(x, 77) • v(x, rj) 



(38) 



ignoring that equation (36) does not fully obey the structure of a 
convolution. Note that we have also restored the prefactor con- 
tained in the weight function. Although being rather symbolic 
equation (38) has strong illustrative power. In the first term we rec- 
ognize the lensing deflection angle a = -2V x O (Bartelmann & 
Schneider 2001). Since the CMB temperature fluctuations are mea- 
sured along the line of sight, only the three-dimensional gradient 
perpendicular to the light ray is considered. Accordingly, the scalar 
product just projects out the velocity components in the plane of the 
sky. Hence, we recover the dipole-like temperature anisotropy pat- 
tern which is usually associated with moving gravitational lenses 
(Birkinshaw & Gull 1983; Gurvits & Mitrofanov 1986). The close 
relation between this effect and the nonlinear iSW-effect has al- 
ready been pointed out by Cooray (2002). However, neglecting the 
divergence part of the velocity field Cooray (2002) does not ac- 
count for the contributions of the gravitational self-energy density 
we discussed before. 

For the computation of the power spectra of the two different 
contributions to the nonlinear iSW-effect we proceed in complete 
analogy to Section 3. In terms of the spectra of the two new fields 
®sg and ®bg the power spectrum reads 



,{k) + tf BGBG (k) + 2fr saBG (k). 



(39) 



The three different spectra involved can be jointly expressed as 



with 



ij>i: 



d/J fxx(k, a,fi)A-(ak) 



fBGBa(k, a, /d) 



A 2 (k^a 2 -2afi+ l) 
a(a 2 - lan + l) 3 ' 2 



2tt 2 - 2an + 1 
a 2 k*(a 2 - 2aLi + l) 2 ' 



/j (2cf 2 + 2a 2 /j 2 - o- 3 yU • 



3afi + l) 



ak^ia 2 — 2afx + l) 2 
and the cross-spectrum 

! (2q-/j 2 ■ 



fsGBaik, a,fi) = - 



■ 3/j + 1 
a 2 ld{a 2 - 2ayL + l) 2 ' 



(40) 



(41) 



(42) 



(43) 



In Figure 5 we plot the three-dimensional power spectra of the grav- 
itational self-energy density and the Birkinshaw-Gull term along 
with the cross spectrum of both effects. In addition, we show the 
resulting power spectrum of the nonlinear iSW-effect and, for com- 
parison, also the corresponding spectrum of the linear iSW-effect 
in the same figure. From scales k > 0.1 (Mpc /hy l on the spectra 
of the © so and SG field, as well as their cross correlation, ex- 
ceed the linear signal by several orders of magnitude. However, 
all three spectra are almost identical, as is demonstrated in the 
lower panel of Figure 5. Only on largest scales the relative differ- 
ence amounts to 10-20%. Therefore, their combination according 
to equation (39) results in the much lower signal of the nonlinear 
iSW-effect, which starts to surpass the spectrum of the linear effect 
on scales k ~ 0.2 (Mpc//?) -1 . Furthermore, we recognize that the 
three spectra peak at k ~ 1 (Mpc//i) _1 , whereas the maximum of 
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Figure 5. Three dimensional power spectra of the gravitational self-energy 
density and Birkinshaw-Gull term constituting the nonlinear iSW-effect as 
well as their cross spectrum (upper panel). Relative strength of the signal in 
comparison to the Birkinshaw-Gull term (lower panel). 




Figure 6. Angular power spectra of the gravitational self-energy density and 
Birkinshaw-Gull term constituting the nonlinear iSW-effect as well as their 
cross spectrum (upper panel). Relative strength of the signal in comparison 
to the Birkinshaw-Gull term (lower panel). 



their combination is shifted to larger scales (k ~ 0.1 (Mpc//z) -1 ). 
Limber projection (cf. equation 28) of the three different spectra 
results in the corresponding angular power spectra shown at the top 
of Figure 6. Since the projection distributes the power of one mode 
over a wide multipole range the differences in the angular spectra 
are even smaller than in case of the three-dimensional power spec- 
tra (see lower panel of Figure 6). From multipoles £ > 1000 all 
three spectra are essentially indistinguishable. 



5 PHYSICAL INTERPRETATION 

To further motivate our interpretation of the contributions to the 
nonlinear iSW-effect we now present some analogies to the theory 
of gravitomagnetic potentials and classical field theory. 

We start with the temperature anisotropy induced by a small 
lens moving with small (transverse) three- velocity v. At lowest or- 
der we have (Birkinshaw & Gull 1983; Gurvits & Mitrofanov 1986) 



! Jdxv 



VO, 



(44) 



where the integral is evaluated along the photon path. Albeit treat- 
ing the involved fields in linear theory, the Birkinshaw-Gull effect 
is obviously nonlinear in the sense that the perturbed quantities en- 
ter as a product. In linear theory the velocity field can be obtained 
as the gradient of the potential (cf. equation 17). Consequently, the 
Birkinshaw-Gull effect is of order (9([VO] 2 ). In perturbation the- 
ory we can naturally construct one more term of this order, namely 
OA®. By virtue of the Poisson equation this is exactly the gravita- 
tional self-energy density we have found in equation (36). 

For simplicity we drop any prefactors in the following. Then, 
both effects can be obtained from the auxiliary vector potential A = 
OVO by taking the divergence 



Ar„, 



2 J d^divA = -2 J A X (OAO + (V<D) 2 ) . (45) 



?CMB 

Rewritting the potential A as 
')VO(x) 



A(x) 



f dV^ 

J-v I* ' 



x'\ 



(46) 



one immediately recognizes its similarity to the gravitomagnetic 
potential 

A GM (x) ee d\x — — - d 3 x ■ (47) 

J'V I* * I J'V I* * I 

The vector potentials can be identified when the change of the 
gravitational potential may be considered as being constant over 
the volume of integration r V. As long as both vector potentials 
coincide, the nonlinear iSW-effect is completely described by the 
Birkinshaw-Gull effect, as well as the gravitational self-energy 
term (at the corresponding order in O) since Schafer & Bartel- 
mann (2006) showed that Ar n i. isw/7cMB = 2 Jd^div Aqm- There 
is an important difference with respect to the theory of electro- 
magnetism. In general, the vector potential describing the magnetic 
field is not fully determined by the scalar potential of the electric 
field. The scalar potential only enters implicitly via the charge den- 
sity. At our level of approximation, however, not only the mass den- 
sity (6 ~ AO) enters the gravitomagnetic potential but the velocity 
field is derived from the scalar potential, too (v ~ VO). Thus, given 
the scalar gravitational potential, the gravitomagnetic vector poten- 
tial is already completely fixed. 

It is interesting to note that the contribution of the self-energy 
density is independent of the assumption VO = const, for x e "~V. 
This can be seen as 



J'V \x-x'\ J v I 



1 



\x-x'\ 



+ ^)V^1 
\x-x'\ 

-4tt6(x)®(x) 



r (b(x') 
+V- d 3 x'd(x')V y ' , (48) 
J-v \x ~ x '\ 



where we have first integrated by parts and subsequently inter- 
changed the order of differentiation and integration. 

The strong similarity to the potential A also indicates why the 
statistics of the Birkinshaw-Gull and the self-energy density term 
are almost identical (cf. Figures 5 and 6). Integrated over a large 
volume one expects the divergence of A to vanish due to Stokes' 
theorem. Accordingly, both terms, considered over a large volume, 
must cancel. Thus, one expects that they reveal the same statistical 
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properties, whereas the small deviations point to the actual differ- 
ences in the two vector potentials A and Aqm- 

Finally, it is worth mentioning an analogy with the theory of 
massless scalar fields. Interpreting the potential O as scalar field 
interacting with an additional scalar field 5 one readily computes 
its energy momentum tensor (cf. Liddle & Lyth 2000) 

7), y = dfl <9 v O - | la,® + 6 oj %v (49) 

where rj^ = diag(-l, 1, 1, 1) and we apply Einstein's summation 
convention. Then the so-called mechanical pressure P of the field 
O is given by one third of the trace of the spatial part of the energy 
momentum tensor, i.e. 

P= \t\= Io 2 -I(VO) 2 -<5<f.. (50) 

We now see that, at our level of approximation, the nonlinear iSW- 
effect measures to some extent the mechanical pressure of the grav- 
itational potential without the contribution arising from its variation 
in time. That the term including time derivatives does not contribute 
is somewhat expected since in our derivation of the nonlinear iSW- 
effect we used the continuity equation to express the time variation 
of the gravitational potential by its variation in space. The analogy 
with the mechanical pressure once more shows that the nonlinear 
iSW-effect is naturally constituted by both the Birkinshaw-Gull ef- 
fect, as well as the gravitational self-energy density. The halo model 
approach of Cooray (2002), however, cannot resolve this close re- 
lationship between these two contributions. 

At the end of this section we shall comment on the appar- 
ent weakness of our approach. Major concerns may arise from the 
fact that we analyse a nonlinear effect starting from linear theory. 
We exploit the commonly used ansatz in perturbation theory that 
a product of first-order fields yields a second-order perturbative, 
i.e. nonlinear quantity. This approach, however, reveals at the same 
time one of the main problems of perturbation theory: it is not clear 
how to generalize this ansatz to higher order. On the other hand, 
this disadvantage is compensated for by the extreme simplification 
provided by this ansatz. It allows for the reduction to one single 
dynamical field, the density contrast or likewise the gravitational 
potential. One can then easily identify the terms which contribute 
at the same order in the perturbative field. Furthermore, the statisti- 
cal properties may be derived most conveniently by tracing back 
higher order correlators to the matter power spectrum (cf. Sec- 
tion 4). Thus, our approach facilitates a fully analytical treatment, 
which in turn offers some direct physical interpretations of the ori- 
gin of the nonlinear iSW-effect. These interpretations are widely 
inspired by analogies and always prefer the illustrative power of 
physical reasoning. 



6 SUMMARY 

In this paper we revisited the theory of the nonlinear iSW-effect. 
We used the continuity equation to express the time evolution of 
the gravitational potential in terms of the divergence of the dark 
matter flux density. The momentum density, being the product of 
density and velocity field, can then be interpreted as source field of 
the nonlinear iSW-effect. 

(i) Our ansatz treats both dynamical fields, density contrast and 
velocity, in linear theory. Consequently, the velocity field is com- 
pletely determined by the density contrast, leaving the latter as the 



only dynamical quantity involved in our analysis. This simplifica- 
tion allows for a fully analytical treatment in contrast to the pre- 
vious work of Cooray (2002) which substantially relies on results 
from the halo model. 

(ii) Exploiting Wick's theorem we derived the angular power 
spectrum of the nonlinear iSW-effect. We confirmed the shape and 
amplitude found by other authors using different approaches (Sel- 
jak 1996; Cooray 2002). Especially, we verified that the nonlinear 
signal surpasses that of the linear iSW-effect at multipoles I > 100. 

(iii) Our analytical ansatz allowed to reveal two contributions 
to the nonlinear iSW-effect: the Birkinshaw-Gull effect (as already 
pointed out by Cooray 2002) and the conformal change of the grav- 
itational self energy density of the cosmic large scale structure. 

(iv) Computing the three-dimensional power spectra and the 
corresponding angular power spectra of the individual contribu- 
tions, we found that they are almost identical, except on large 
scales, and their individual detection is certainly impossible for on- 
going CMB experiments. 

(v) We showed by simple arguments from perturbation theory 
that it is natural that the Birkinshw-Gull term is accompanied by 
the gravitational self-energy density term. Furthermore, we pointed 
out several analogies with the theory of gravitomagnetic potentials 
and scalar fields to highlight the physical meaning of the different 
contributions to the nonlinear iSW-effect. 

Despite its confinement to linear theory our ansatz allows for a 
deeper understanding of the physical processes underneath the non- 
linear iSW-effect, especially regarding its relation to other phenom- 
ena like the Birkinshaw-Gull effect. 
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